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Abstract 

We consider the dynamics on a quantum graph as the limit of 
the dynamics generated by a one-particle Hamiltonian in M. 2 with a 
potential having a deep strict minimum on the graph, when the width 
of the well shrinks to zero. For a generic graph we prove convergence 
outside the vertices to the free dynamics on the edges. For a simple 
model of a graph with two edges and one vertex, we prove convergence 
of the dynamics to the one generated by the Laplacian with Dirichlet 
boundary conditions in the vertex. 

I. INTRODUCTION 

Several physical systems, like electronic nanostructures or periodic solids 
made up of aromatic molecules, display the common feature that the motion 
of the electrons can be thought of as being confined in one or more directions 
by a strong potential barrier, which prevents the particles from escaping from 
the structure and allows propagation in the remaining (free) directions. 

In solid state physics, a concrete example of this phenomenon is given 
by quantum wires and carbon nanotubes, where the high purity achieved in 
fabrication techniques and the weakness of electron-phonon interaction give 
rise to ballistic transport. Therefore, using the strong-coupling method, one 
can in first approximation model the interaction of electrons with the crystal 
assuming that they move freely with an effective mass m (a more detailed 
discussion of the physical hypotheses on which this approximation is based 
can be found in Duclos and Exner (1995), Londergan et al. (1999)). 
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Taking into account the strong potential barrier which keeps them con- 
fined, their dynamics is then given by the one-particle Hamiltonian 



where we have chosen suitable units so that h and the effective mass are 
equal to 1. 

The parameter e is the natural small parameter of the problem and is 
linked to the ratio l/L, where L is the characteristic length of the wire along 
the free direction where the electrons can propagate and I is the analogous 
length in the confined directions. 

In general, the function W is assumed to be zero on the quantum wire 
and strictly positive outside, so that when e becomes small one expects the 
electron to be better and better confined to the wire. For this reason, W is 
called the constraining potential. 

Our aim in this paper is to analyze the dynamics generated by (£[J in 
the limit e — > for two-dimensional systems constrained to a singular one- 
dimensional manifold given by a graph F. 

We consider an explicit form of the constraining potential, i. e., the dis- 
tance of a point from the graph. This potential is obviously continuous, 
but does not belong to C 1 (M 2 ), because of the vertices. Outside the ver- 
tices, the potential is a quadratic function of the coordinates, so it can be 
regarded as the first non zero term of the Taylor expansion of more general 
functions which are zero on the graph, and whose gradient is zero on the 
edges. This is the natural class of constraining potentials considered in the 
literature to model holonomic constraints in classical and quantum mechan- 
ics (Bornemann (1998), Froese and Herbst (2001), and references therein). 
We will argue that the results we find can be generalized to a wider class 
of potentials, whose Taylor expansion near the graph contains higher order 
terms. 

As it is clear from the brief remarks given above, this case is relevant 
for solid state physics in order to determine the leading order behavior for 
the dynamics of electrons moving in a branched nanostructure, made up of 
several wires which meet in a crossing region represented by the vertices of 
the graph in the limit e — > 0. 

Another important field where equation (0) can be applied is given by 
theoretical chemistry, in the so-called Quantum Network Model (QNM) (for 
a recent review which contains also some comparison with experimental data 
see Amovilli et al. (2004) and references therein). This model is used to 
study the motion of valence electrons (also called 7i~-electrons) in aromatic 
molecules or periodic solids like graphene. In first approximation, they are 
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Figure 1: Skeleton of naphthalene molecule 



thought to move freely through the skeleton of the molecule determined by o~- 
electrons, which create a potential keeping 7r-electrons confined to the molec- 
ular structure. The most famous example, described in figure [TJ is probably 
naphthalene molecule, first studied by Ruedenberg and Scherr (1953). 

The characteristic feature shared by the examples presented is that we 
always expect the limiting dynamics to take place on a graph, on which we 
have to define a suitable Hamiltonian. The couple made up of the graph 
and a differential or pseudo-differential self-adjoint operator defined on it is 
usually called quantum graph (a detailed review on this topic is given in 
Kuchment (2002, 2004, 2005)). 

The main problem one runs into if one wants to describe the physical 
systems mentioned above by a quantum graph is that there are a host of 
self-adjoint Hamiltonians which can be defined on it. 

If, for example, one assumes that the dynamics outside the vertices is free, 
i. e., that the Hamiltonian is the Laplacian —d 2 /dx 2 , where x is the natural 
arc length coordinate on the edges of the graph, and that the graph has 1 
vertex and n edges, then the possible self-adjoint extensions are determined 
by n 2 parameters (Kostrykin and Schrader (1999)). 

To eliminate this ambiguity, a natural procedure is to consider the quan- 
tum graph as a limit of a more realistic model. We call the method employing 
a strong constraining potential a soft approximation of the graph. As we have 
pointed out above, this is a physically reasonable approximation, which, for 
a smooth constraint, has been investigated, e. g., in Belov et al. (2004), 
Dell' Antonio and Tenuta (2004), Froese and Herbst (2001). As far as au- 
thors' knowledge is concerned, the remark by Kuchment (2002) that "the 
graph case (i. e., in the presence of vertices) has not been explored" is still 
largely true. 

Another appropriate choice, which we call rigid approximation of quan- 
tum graphs, is to consider a "thickened graph", composed in the simplest 
cases of thin tubes of radius e, which has the same topology as the original 
graph and reduces to it when e goes to zero. It is reasonable to suppose that 
the motion of the electron in this thickened structure is free, but one needs 
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to specify boundary conditions to get a well-defined Laplacian. The most 
natural choice is to use Dirichlet boundary conditions, which correspond to 
an infinite constraining potential barrier. Some light on this case has been 
shed only very recently by Post (2005), who considered thickened graphs 
which are strictly smaller than the ones defined by the distance function, 

V £ , r :={(x,y)eR 2 :d r (x,y)<e}. (2) 

This is another reason why the study of a constraining potential given by 
the distance is interesting. 

For technical reasons, much more attention has been devoted to the case of 
Neumann boundary conditions, which is by now well understood (Exner and 
Post (2005), Kuchment and Zeng (2001, 2003), Rubinstein and Schatzman 
(2001), Post (2006)). 

With the exception of Post (2006), all the papers we have quoted deal 
with the convergence of the spectrum of the Laplacian, defined on a thickened 
graph, to the spectrum of a self-adjoint operator, defined on a graph whose 
edges have finite length, while we are primarily interested in the convergence 
of the dynamics and eventually to scattering theory. 

A weak form of resolvent convergence for Neumann rigid approximations 
was studied by Saito (Saito (2000, 2001)), but his results do not allow to 
infer the structure of time evolution. 

In Post (2006) these results are improved, and norm resolvent conver- 
gence is established, under the hypotheses that the vertex neighbourhoods 
are "small" (for the precise meaning of the term we refer to the original pa- 
per) and, for graphs embedded in R 2 , that the angle between two different 
edges has a global lower bound. 

Another important difference is related to the class of initial conditions 
we consider. From a physical point of view, one expects that, when e — > 0, 
the Hamiltonian gives rise to fast oscillations of the electron in the direc- 
tions orthogonal to the edges of the graph. To prove the results mentioned 
above, one projects, roughly speaking, on the ground state of this trans- 
verse oscillation. Nonetheless, in the spirit of adiabatic perturbation theory 
(Teufel (2003)) we expect to be able to get an effective dynamics inside every 
transverse subspace, because they become broadly separated in energy when 
e goes to zero. For this reason, we consider initial wave functions which 
are localized inside one edge, and belong to an eigenspace of the transverse 
Hamiltonian (which will be defined more precisely in section II). 

We consider longitudinal initial conditions which are independent of e. 
This corresponds to study longitudinal states which vary over a wavelength 
which is much bigger than the transverse one. As it has been stressed in 
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Belov et al. (2004) however, longitudinal states are not homogeneous, and it 
would be interesting to consider also wavefunctions which vary on a scale of 
order e 1 ^ 2 for example, analyzing in this way a semiclassical limit. 

One should also remark that, in the Neumann case, the energy of the 
transverse ground state is independent of e, because the Neumann Laplacian 
has always the eigenvalue zero corresponding to the constant function, while 
in the Dirichlet and soft approximation cases, the energy of every transverse 
mode tends to infinity when e — > 0. This makes it necessary to subtract a 
divergent phase to get a finite result. 

We will now describe briefly the structure of this paper. 

In section II we study the convergence of the unitary group generated 
by (0) for an arbitrary graph (i. e., with an arbitrary number of vertices 
and edges). As we have already mentioned, we consider initial conditions 
which belong to a transverse eigenspace and are localized inside one edge of 
the graph and we choose a constraining potential given by the distance of a 
point from the graph. Using weak convergence methods we study the limit 
flow on the graph outside the vertices. To describe completely the limit flow, 
we must study its structure in a neighborhood of the vertices. 

This may be a difficult task; this can be seen from the exact treatment 
we give in section III of a simpler system, in which the graph is a continuous 
curve in the plane. 

In this example, the graph has one vertex and two straight edges at an 
angle < $ < ir. In this case we approximate the graph by a sequence 
of smooth curves converging to the graph when e — ■> 0. We prove that 
generically (in particular if the curvature of the approximating curves is ev- 
erywhere non-negative) the limit dynamics along the graph correspond to 
Dirichlet boundary conditions at the vertex. 

The proof is achieved by reducing the problem to the study of the dy- 
namics with Dirichlet boundary conditions on the boundary of narrow tubes 
containing the graph, using a refined version of a theorem due to Froese 
and Herbst (2001). This however is not possible for every smooth curve ap- 
proximating the graph, and a condition on the curvature comes in. From 
this result, it seems that the constraining potential and Dirichlet boundary 
conditions are not always interchangeable, as one could naively think. 

Even though the geometry of the graph is very simple, this model demon- 
strates a mechanism through which adiabatic decoupling among different 
transverse modes takes place. In particular, it shows that the bound states 
localized near the vertices that can arise (and indeed do arise if instead of the 
constraining potential one considers a narrow tube with Dirichlet boundary 
conditions, see Carini et a/.(1992, 1993), Goldstone and Jaffe (1992)) do not 
interfere with the propagation of product states localized inside one of the 
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edges at the leading order, because their spectral distance becomes infinite 
in the limit. 



II. CONVERGENCE OUTSIDE THE VERTICES 

In this section we consider a finite metric graph, denoted by T, whose 
edges can have infinite length. We assume that it is embedded in R 2 and, 
for the sake of simplicity, that all the edges are straight lines. 

We denote by V = {fj}i e / the (finite) set of vertices and by E = {ej}j e j 
the (finite) set of edges connecting them. We assume that there are no 
isolated vertices. 

A graph is said to be a metric graph if to each edge e is assigned a 
length l e e (0, +00]. Edges of infinite length arise naturally if one considers 
scattering theory on graphs (see, e. g., Melnikov and Pavlov (1995) and 
references therein). 

We can now identify each edge with a finite or infinite interval [0, l e ], with 
the natural coordinate x e along it. One can also define function spaces (e. 
g., L p spaces, Sobolev spaces); in the case of Sobolev spaces, one must have 
some care at the vertices (see Kuchment (2004)). 

As mentioned in the introduction, we approximate the dynamics on the 
graph using an Hamiltonian, acting on L 2 (R 2 ), with a constraining potential 
given by the square of the distance from T: 



One could use a different potential, whose Taylor expansion away from the 
vertices contains higher order terms, but we make the important assumption 
that the Hessian is constant along the edges with the same value on all 
edges. This condition is reasonable in view of the analysis in Dell'Antonio and 
Tenuta (2004) and Froese and Herbst (2001) when the system is constrained 
to smooth submanifolds. 

In this case, the energy of the transverse oscillation appears as a poten- 
tial energy in the longitudinal motion, in the form u(x)/e, where u is the 
frequency of the oscillation. If u does not depend on x, the resulting phase 
factor in the dynamics can be discarded. Otherwise, it originates a constrain- 
ing potential along the edge, so that in the limit e — > we expect that the 
wave function concentrates along the minima of u. 

We denote by U t (e) the unitary evolution associated to Hamiltonian (jSJ), 




q,qe 
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(3) 



U t (e) := exp(-itH(e)), 



(4) 
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and we take as initial a state which "lies in a subband", i. e., is in a fixed 
transverse mode and localized within one edge. These are the states which are 
thought to describe the propagation of particles in semiconductor structures. 
We have then 

(x ej0 , y ej0 ) = / (x ejQ ) K (y ej0 ) , (5) 

where / G CJf (0, l e . ) and $^ is an eigenstate of the harmonic oscillator, 



1 d 2 



1 



2dy 2 + 2e 2V 
1 

= n+-. 



(6) 



C J0 



is the natural coordinate along the edge ej and y ej is the corresponding 
coordinate in the orthogonal direction. x ejQ varies in the interval [0,/ ej J (or 
[0, +oo) if the edge has infinite length) and, since the edges are straight lines, 
y ejo is well defined and assumes values between — oo and +oo (to simplify 
the notation, from now we denote these coordinates just by Xj and yj). 



Figure 2: Schematic representation of the initial state 




Applying U t (e) to ip we expect the appearance of a strongly oscillating 
factor, given by exp(— iE n t/e). To avoid this (irrelevant) phase, we consider 
the modified unitary group: 

H(e): Z H(e) - f , _ (?) 
ift := U(e)ip := exp(-itH(e))tp - 

To analyze the adiabatic decoupling, we split ty\ according to the different 
transverse components in each edge: 



8?(t, Xj ;e):= J dy 3 $ £ m ( % )*^(a;„ y,) 1 , (8) 

where Xj is the natural coordinate along ej and yj is orthogonal to it. 

1 With an abuse of notation, we denote by ipt( x ji Vj) the function ^| written in coordi- 



7 



Proposition 1. Let P™ the operator from S(M?) to H k (R) defined by 

P^{x 3 ):= I dyj&MWxjM), (9) 

then Pp extends to a unique operator (of norm 1) from H k (M. 2 ) to H h 
for every k G N 7 k > 0. 

Proof. Given G S(R 2 ) it is clear that 

d l x (P^)(x) = [ dy^ m (yrd l Mx,y)=<^ m ,d l x ^> L 2 {Ry) 
Jr 

=► \d l x (P^)(x)\ 2 < \\* m \\h<R) ■ J dy\d l J(x,y)\ 2 
=► J dx\dl(Pp^)(x)\ 2 < ||<f m ||! 2(R) -||^||i 2(R2) 



wiiW> = En^ m ^nW) ^ iim 2 £ii^(*,v)iiw ( 10 ) 

□ 



Z=0 1=0 



Corollary 1. The components s™(t,Xj;e) are well defined, belong to H 1 
in the variable Xj and satisfy 

sup \\s^(t, ■;e)\\ L 2 {R) < const. . (11) 
te[o,r] 

Proof The domain of the quadratic form associated to H(e) is given by 

Q(H(e)):=H\R 2 )nQ(d 2 ) , (12) 



where Q(d 2 ) := {ip G L 2 (W) : d r (x,y)ip G L 2 (R 2 )}. ^ belongs to Q(H(e)), 
so V* is in ^(R 2 ). □ 

Lemma 1. 

||-^(^)V'o|| < C (independent of e). (13) 



nates (xj ,Uj). Since the different systems of coordinates associated to each edge are linked 
to one another by a rigid motion of the plane, this does not modify the differentiability 
or integrability properties of ipf. Note, instead, that is an eigenfunction of the 

harmonic oscillator in the variable yj. 
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Proof. Since the Laplacian is invariant by rotations and translations, we have 
(for simplicity we drop the index j in x and y) 

H(e)^ = -~d 2 J •<!>„--/• d^ n + ^4f ■ $ n - ■ $ n = 

2e 2 e 
= -\d 2 J-<S> n + ±- 2 (4-y 2 )f.<S> n . 

Let a be the infimum of the support of /, and b > a the supremum. 
Since / is supported inside the edge e J0 and near each edge the distance from 
the graph is equal to \yj\, the function ^2 (4 — y 2 )f • $n will be zero when 
\y\ < D, where D, depending on the support of /, can be small, but it is 
strictly positive. Therefore we have 



f dx I dy -y 2 ) 2 \f(x)K(y)\ 2 - (14) 



Now we use the following two properties: 

• <ip is equal to a polynomial of second order in the variables (x,y); 

• $ £ is equal to a polynomial in y/e 1 ^ 2 times exp(— y 2 /2e). 

The norm (j!4j) contains then terms of the form (P and Q are polynomials) 



dx \f(x)\' 



dy P(x, y)Q(x, y/e 1/2 ) exp 



\v\>d 



y_ 

e 



dx 



dy P(x,y)Q(x,y/e 1/2 )exp 



\y\>D 



M)y 



exp(-My 2 ) < 



< exp 
= 0(e~ c / £ ). 
This implies that 



dx 



dy P(x, y)Q{x, y/e 1/2 ) exp(-My 



y\>D 



\^{d 2 T -y 2 )f^ n \\ 2 = 0{e-^) 



and therefore the thesis is proved. 



(15) 
□ 



9 



□ 



Corollary 2. For every system of coordinates (xj,yj) associated to an edge 
Cj we have 

l\\d Xj m 2 + \WdyM + ^IK^II 2 < j- (16) 

Corollary 3. Let F dr >$ be the characteristic function of the set {{x,y) : 
dr{x,y) > 5}, where 5 is any positive number, then 

sup \\F dr > s ^ £ t \\m^) = 0(e 1/2 ). (17) 
te[o,T] 

Proof. Last term in ()16|) gives 

\\\drr t \\ 2 <Ce. 

Therefore we get 

6 2 ■ < i/4, F dr > s i> £ t >< < dr i>l F dr > 5 d r ^ t >< 2Ce. 

We are now ready to prove 
Theorem 1. For m ^ n and j G J 

sf(t, Xj ;£) A 0, e^O, (18) 
where the convergence is in the weak* topology of L°°((0, T), L 2 (IR)). 

Proof. For convenience of the reader, we recall that a bounded sequence f e 
of functions in L°°((0, T), L 2 (R)) converges to a limit /„ G L°°((0, T), L 2 (M)) 
in the weak* topology if and only if 

dt lf(t) < X('), fe(t, •) > L 2(R)^ / dt tp(t) < *(•), fo(t, •) >L2(R) 

JO 

for every function ip G L x (0,T) and every function x £ L 2 (R). 

It is a standard fact about weak* topology that it is enough to consider 
only if and x hi dense subsets of I/ 1 ((0, T)) and L 2 (R) respectively (see, e. g., 
Bornemann (1998), Appendix B or Rudin (1973), theorems 3.15 and 3.16). 

Let us then show first that 

sup \<x(-),sf(t,-;e)> L 2 (R) \^0, e -> 0, (19) 
te[o,T] 
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for every function x £ C^(IR\{0, L.}) if the edge has finite length and for 
every x e C[5°(IR\{0}) if ej has infinite length. We consider explicitly only 
the latter case, the former being analogous. 

If the support of x is contained in (— oo, 0) then 

l< X,sf >L\R)\ = \<X-$ £ m ^l >L2(R2)| = \< X -& m ,F x $\ >| 

<||x-^M| j p x ^II = c'(£ 1/2 ) ) 

where F x is the characteristic function of the support of x an d we have used 
corollary El 

If the support of x i n contained in (0, +00), then, following the proof of 
lemma we can show that 

H(£)[x(xjW m (yj)} = (we drop the index j) - ]-d 2 x x • ® £ m + 

Since x is supported away from the vertex (located at Xj = 0), an equation 
similar to (|15j) holds: 

11^(4 -y 2 )x-® £ J 2 = o(e-^). 

Since m 7^ n, we have then 
This implies 

< x(0.*™(v;<0 >=< xi-Wm^t > = ^— < xi-WmMzWt > + 

J m — n 

+ < \dlx ■ ^ # > +0(e^) - 

m — n 2 

= (lemma GJ 0(e). 

Now, if cp G L 1 ((0, T)), we get 



< y\\ L i- sup 1 < x(-)> s f (*.-;e) > 

te[o,T] 



but we have just shown that the right-hand side goes to zero for x G C^°(IR\{0}) 
(or Q5°(IR\{0, Z e }) for an edge of finite length) which is dense in L 2 (R). □ 
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Theorem^ shows that, although in a weak sense, there is indeed adiabatic 
separation between the different transverse states even in the presence of 
vertices, if the initial state is localized in two senses: first, it has to be 
localized inside one edge to avoid mixing between the different transverse 
states associated to each edge and second, it has to be in one (or a finite 
number of) transverse band(s). 

Since the limit of sj 1 for m 7^ n is zero, to analyze in a complete way the 
(limit) evolution of ipo we have to determine the behaviour of s™, j G J as a 
function of time. 

Theorem 2. There exists a weak* convergent subsequence of s^(t,Xj',e) in 
L°°((0, T), L 2 (M)) (denoted again by the same symbol), whose limit s"(i, Xjj 0) G 
L°°((0,T),L 2 (M)) satisfies 

id t s?(t,x j] 0) = -±%s](t,x j ',0) m £>'((0,T)x(0,Z e .)). (20) 

Remark 1. By corollary^ Xj '. <s) is db bounded sequence in L°°((0, T), L 2 (IR)). 
Since the balls in L°°((0, T), L 2 (R)) are compact metric spaces with respect 
to the weak* topology (see the theorems in the book of Rudin quoted above), 
a weak* convergent subsequence certainly exists. 

Moreover, if one shows that all the weak* convergent subsequences con- 
verge to the same limit, then this implies that the sequence itself converges. 

The equation satisfied by the limit in theorem 121 is clearly independent of 
the subsequence, but it does not determine the behaviour of the limit in the 
vertices, so we cannot conclude convergence of the sequence. 

For this it would be necessary to control the behavior of the sequence 
in a neighborhood of the vertices. This difficulty (which is not present for 
smooth submanifolds, where the same strategy has been successfully applied 
by Bornemann (1998) in the classical case) is linked with the fact that the 
operator —d 2 /dx 2 defined for functions which vanish in a neighborhood of 
the origin has many self-adjoint extensions which define different dynamics. 

We split the proof of the theorem into a number of lemmas. 

Lemma 2. s"(£, Xj',e) belongs to C 1 ([0, T], L 2 (M)) and moreover it is an 
equicontinuous sequence of function from [0,T] to L 2 (R). 

Proof. Let us denote by s](t,Xj-,e) G C°([0, T], L 2 {R)) the function 
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Using proposition ^ we have that 



s?(t + h,-,e)-s*(t,-;e) 



< 



h 

U t+h (e) - U t (e) 
h 



-S?(V;e) 





pn 

3 


L 2 (R) 





U t+h (e) - U t (e) 



+ iU t (e)H(e) 



+ iU t {e)H{e) 



o. 



This proves that 

id t s](t,x j; e) = s](t, Xj ;e) = P™[-iU t {s)H{e)^\. (21) 

Since (s^oll is bounded ( lemma QJ), 1 1 s J 2, (t , xj] e)\\ is bounded, therefore 



\sm,-;e)-s?(t',-,e)\\ LH 



drd T s™(T,-,e) 



< 



dr iia^r,-;*)!! <C\t-t'\, 



showing that s^it.x^e) is an equicontinuous sequence. 

Corollary 4. There exists a subsequence s™(t,Xj;£) which satisfies: 



□ 



1. sY{t,Xj;e) converges, in the weak topology of L 2 (R), uniformly in t, to 
a limit s™(t, Xj] 0) G L 2 (R). Moreover, the limit is continuous in t in 
the weak topology of L 2 . 

2. d t sf{t, Xj ;e) A d t sf{t, Xj ; 0) in L°°{ (0, T), L 2 (R)), where the deriva- 



tive df-S^it, Xj] 0) is to be interpreted as derivative in P'((0, T) 



x 



Proof. The sequence sJ l (t,Xj;e) is contained in a ball in L 2 (R). This ball 
is a compact metric space with respect to the weak topology. Since the 
sequence is equicontinuous with respect to the strong topology, it will be 
equicontinuous with respect to the weak topology too. Therefore, the theo- 
rem of Ascoli-Arzela (see, e. g., Royden (1988), theorem 10.40) proves the 
first point. 

Equation (J2T|) implies that d t s™{t, Xj] e) is a bounded sequence in L°°( (0, T), L 2 (M)), 
so, extracting possibly another subsequence, we have that there exists gj 1 e 
L°°( (0, T), L 2 (M)) such that (again, we denote the subsequence with the 
same symbol as the sequence itself) 



d t s^{t,Xj;e) A g 



m 
3 ' 
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but this implies that V (p e Cg°((0,T)), V \ e C, 

f -T 



OC 1 



r 



dt / dx gf(t,x)ip(t) X (x) = 
Jr 

dt<p{t)< X (-),9?(t,-) >l*(r)^- / rft^(t)<x,^ m (t,-;e)>= 

/ dtcp(t)d t <x,sf{t,-;e)>= 
Jo 

-[ dtd t <p<x,s?(t,-,e)>^- f dtd t <p< X ,s?(t,-;0)>= 
Jo Jo 

dt / dx sf(t, x;0)d t <p(t)x(x), 
Jr 

gf = d t sf(t,xj;0) in V'((0,T) x R x ). 

□ 

We can now prove theorem |21 

Proof. We suppose that edge ej has infinite length. The proof for an edge of 
finite length is analogous. 

Corollary El together with the proof of the first part of the proof of the- 
orem ^implies that 

sup \< X ,sf(t r ;e)>\ = 0(e 1 / 2 ), 
te[o,T] 

for all j G J and for all x £ C§° ( — oo,0), but the first point of corollary |U 
gives 

< x, s?(t, •; 0) > L 2 (M) = lim < x, sf(t, •; e) >= 0. 
Equation (J2TJ) allows us to write, for all x e C^°(0, +oo), 

< X, ^ t s™(i, •; e) > L 2 (R) =< x • K, H{eW t > L 2 (R2) =< £T (e)x • K, i>\ >= 

=< ~\dlx ■ K, # > + < - Z/ 2 )X ■ ^ £ >=< ~2 d lx, s?(t, ■; e) > L2(M) + 

+ <^(4-y 2 )X^n>^> 

Since x is supported in (0,+oo), equation (|THj) holds also in this case, 
therefore 

11^(4 -y 2 )x-K\\ 2 = 0(e-^). 
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We have then, for all ip e Cg°(0,T), and for all X e Cg°(0, +00), 

/ dt <p(t) < x, id t s](t, •; 0)) > L 2 (R) ^ [ dt ip(t) < X , id t s](t, ■; e) > L 2 (R) = 
Jo Jo 

= J dt <p(t) < -^d 2 xX ,s](t,-;e) > L2(R) +O(e- c / £ )y\\ L i {0 , T) 
- J dtip(t)<-^d 2 xX ,s](t,-,0)> L 2 {m) , 

td t s](t,x 3 ;0) = -±%s](t, Xj] 0) in V'((0,T) x (0,+oo)). (22) 

□ 

III. A GRAPH WITH TWO EDGES 

In this section, we are going to put forward a different kind of soft ap- 
proximation for a graph with one vertex and two infinite edges. We denote 
it by r#, where d is the angle made by the two edges, < d < ir. 

As we have already said in the introduction, we do not consider directly 
r,?, but we approximate it by smooth curves, r^, whose curvature becomes 
bigger and bigger in a region whose width, given by 5, goes to zero and we 
consider a potential constraining to this family of curves. 

More precisely, to specify the approximating curves we need only to spec- 
ify their curvature, kg, because, as it is well known, this determines the curve 
up to rigid motions of the plane. Naturally, we want that, when 8 goes to 
zero, the curves tend to the graph. This in particular implies that the turning 
angle has to become equal to 1? when 5^0. 

A simple choice which satisfies these requests is (s is the arc length pa- 
rameter) 

k s (s) := ^ k (^) » J ds k ( s ) = 1» 

(23) 



fceCS°(-l,l), 



k = l \s\ < 1/2 
k = Q |s|>3/4, 



which amounts to deformate the graph in a neighbourhood of the vertex 
replacing it with an arc of a circle. Note that the 5 scaling is fixed by the 
request that the turning angle of the approximating curves be 
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Actually, our result does not depend on this specific choice we have made, 
because, from the proof, one can see that the only essential ingredient is the 
singularity 1/5, which is forced by the requirement that the turning angle is 

We consider the Hamiltonian 

H(e, 5(e)) = ~A + \w s(e) , 5(e) -> when e -> 0, 
where, for simplicity, we suppose that 

W s{e) (x,y) = -d\ e) (x,y) = -dist[(x, y), F^e)] 2 - 

The remark we made above about the possibility to generalize the analysis 
to potentials with constant Hessian, applies here too. As in the previous 
section, we are interested in the time evolution of a product state which is 
initially localized away from the vertex. 

We expect that the particle oscillates very fast along the direction nor- 
mal to the curve, so to analyze the motion we should use a suitable system 
of coordinates, adapted to the curve. A natural choice is given by tubu- 
lar coordinates, which are a set of local coordinates suited to study tubular 
neighbourhoods of embedded submanifolds (see, e. g., Lang (1995)). If the 
submanifold one considers has codimension (and dimension) bigger than one, 
then the metric in tubular coordinates is in general not diagonal, and the 
Laplacian in these coordinates contains a gauge term which couples the lon- 
gitudinal and the transverse motion (Mitchell (2001) and references therein). 

In our case, since both the dimension and the codimension are equal to 
one, these problems do not appear. 

Given a smooth curve C, with parametric equation ( : Q — > M 2 such that 
|9 s C(s)| = 1, we can describe the position of points in a tubular neighbour- 
hood N of C via the curvilinear coordinates (s, u) defined by 

q(s,u)=C(s)+un(s), (24) 

where q is an arbitrary point of N, n(s) is the normal unit vector to the 
curve and u is assumed to be smaller than the radius of curvature. 

In the case we are dealing with, this means that curvilinear coordinates 
are defined only in the region 

{(x,y)eR 2 :d s (x,y)<Q}, (25) 

where g is the radius of curvature of r^ j( j. When k$ is different from zero, this 
quantity is proportional to 5 itself, so, by hypothesis, it goes to zero when 
e -> 0. 
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To get rid of the region {(x,y) : dg(x,y) > 5} we will use a theorem, 
proved first by Froese and Herbst in the more general context of a potential 
constraining to a submanifold (proposition 8.1 in Froese and Herbst (2001)), 
which basically says that if one starts from an initial state more and more 
localized near the constraint, then all that matters for the time evolution is 
a small region near the constraint itself. We repeat the proof of Froese and 
Herbst because we need to keep track of the dependence of all the constants in 
the estimates on 5, to apply them to the region {(x, y) G M. 2 : ds(x, y) < 5(e)}. 

Theorem 3. Let if) G L 2 (R 2 ), \\ip\\ = 1 and \\H(e,5)ip\\ < (C x indepen- 
dent of 5). Then, 



||iV 5 e-^*V|| < (2C 1 ) 1 / 2 ^. (26) 

o 

F. indicates the characteristic function of the region indicated. 

Moreover, let Ho(e, 5) be the Hamiltonian H(e, 5) with Dirichlet boundary 
conditions on the set {(x, y) G M 2 : dg(x,y) = 5}. 

Let us suppose that 5 = 5(e), \im e _, 5(e) = 0. Taking into account (|26[). 
let us also assume that 

lim — = 0. (27) 
5(e) 

Then, for all t G [0, T], we have 

- e- itS °^F ds<s iP\\ < C 2 (C U T) ^ + ^ . (28) 
Remark 2. The theorem implies that if we choose a 5(e) such that 

1/10 

lim — - = (29) 
5(e) 
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then we can restrict ourselves to analyze the Dirichlet Hamiltonian H^(e, S(e)), 
which is localized inside the region where tubular coordinates are defined. 
This, however, means that we have to consider a "tube" encircling the graph 
whose diameter is much bigger than the localization of the transverse states, 
which for an harmonic oscillator is e 1 ^ 2 . 

Remark 3. As already observed in Froese and Herbst (2001), the estimate 
is not optimal. 



Proof. Let us first prove (|2*H|) . 

Since S)ij)\\ < we have from Schwarz inequality 

< ip,H(e,5)ip >< — . 
This implies immediately that 

>< >= 1|| VV ,||2 + 1 ||^||2 ^ < 2Ci£ _ m 

e I le A 

It follows then 

5 2 < F ds >si),F ds > 5 i> > < < d 5 F ds >s^,d s F ds >s^ > < \M\\ 2 < 2C l£ . (31) 

The same argument can applied also to e~ %tH ^ £ ' 5 'ij), so (j2T?j) is proved. 

We need now to prove an estimate on the behaviour of the gradient of i/j 
away from the graph. 

Let x e Cg°(M) be 1 when 1/4 < |ar| < 3/4 and when |x| < 1/8 or 
|x| > 7/8, then the function 

( \ ~( 1 i u 1 A i 



2(a-Ai)' 1 4 2(a-Ai) 

will be 1 when Ai < \u\ < a and for \u\ near zero. If we choose Ai and a 
such that < Ai < a < 5, then x is we U defined (and G C^°(R)) when u is 
the coordinate along the direction normal to the curve T^j. 
We have then 



a Vii\\ = \\F Xl<d5<a V( X m<\Wxi>)l 



Using again the Schwarz inequality and the fact that x'P ^(A) (the po- 
tential is bounded on the support of x) we get 

iiv(^)ii<iiA( X ^)ir /2 iix^ir /2 , 
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so, to estimate ||i 7 Ai<d 4 <aV , ?/>|| we need to get an estimate on ||A(xV')||- To 
obtain it, we use an energy estimate of second order, i. e., we calculate the 
quadratic form associated to H(e,5) 2 . 

H(s,Sr = \\p\* + (^) 2 + E^2?fe-^ A ^ ( 32 ) 

where p = — iV. The first three terms are positive operators, while if we take 
the mean value of the last one with respect to the state x4> we get 

< X?P,Ad 2 s x4> >= / dxdy \x^\ 2 Ad 2 = / dsdu [1 + uk s (s)] \x4>\ 2 [1 + ukgis)]' 1 - 

Jd s <S J\u\<6 

■d u {[l + uk 5 (s)]2u} = / dsdu[l + uk 5 (s)}\x^\ 2 -2+ dsdu\x^\ 2 2uk s {s) < 

J\u\<8 J\u\<5 

<C\\x^\\ 2 , 

where the Jacobian of the change to curvilinear coordinates is given by 

and in the last step we have used the fact that 

sup \uks(s)\ < 5k$(s) = $k( - ) < const, (independent of 5) 



\u\<5 \° J 

Taking the mean value of (|32|) with respect to x^P we obtain then 

\\U( X n\ 2 <\\H(s,5)(xm 2 + ^ 

which can be written equivalently as 



^||A( X ^||<^ + ||%,^|| + ^||[A, 
The last term is equal to 

[A,x]^ = (A X )^ + Vx- V^, 
and we can estimate its norm changing to curvilinear coordinates, 

1 u 

Vx = d x x{x{u))— — — n(s), 

2{a — Ax) \u\ 

A X = (1 + k^)-%[{l + k s u)d uX ] = d 2 uX + (1 + k s u)- x k s d uX 
= dlx(x(u))— — ?— ^ + d x x(x(u)) 



2(a-Ai) 2 X/Xy y /J 2{a- \ x ) \u\ 1 + kgu' 
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Using (J3*U|) to estimate HV^H, we have then 

II [A, x\M\ < cl/2( ° l , + ( C , n 2 + 77 TT - (34) 
e A '^(a — AiJ (a — Ai) z o(a — AiJ 

In what follows, we will need to choose a and Ai proportional to 5. Assump- 
tion (|27[) implies then that all terms in (|34|) are at most of order e~ l . 
To sum up, we have 

I|A(X0)II<§, (35) 

from which it follows that (assuming that a; and Ai are proportional to S and 
that 5(e) satisfies flZ3)) 

||*W«W|| < Cs-^ 2 = -^j- 2 . (36) 

Let now £ be a function in C^(R) such that £(x) = 1 when \x\ < 1/4 
and £(x) = when \x\ > 1/2. We define the function £ by the equation 
£(«) := £(u/5), where u is the curvilinear coordinate normal to the curve. 

Because of (|2*^jl . to prove (|2*5|l is enough to show that 



/c-l/4 F l/2 

H^^^-ftH^^ _ ev ,|| < c 2 (C 1 ,T) + £ — 
for i G [0,T]. Let 

: = e <tHu(e,i)^ e -ttH( e ,«)^ _ 
Integrating the derivative we have 

<t>t,e,s = i [ ds e %sn ^ 5 \H D (e,5)i~ iH(e,5)]e- lsfl ^ s ^ = 
Jo 

= [ ds e isSD{£ > 5) [V£ • p - (i/2)A£}e- isS(£ > 5) ij, 
Jo 

therefore 

Ht,s,s\\ 2 = f ds <e-^^0^,[Ve-p-(V2)A£]e-^V>. 
Jo 

Let now ( be a C^°(R) function which is 1 on the support of d x £ and when 
\x\ is near zero. As above, we denote by ((u) := ((u/5). We can then write 

II^Mll 2 < f ds ||Ce- iS ^ M) ^ll(l|V£-pe-* 5 )^|| + ||(l/2)A£e-^ 5 )^||) < 
Jo 



(37) 
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where we have used (I36|) and the definition of £. 
Now 

< fa^g, H D (e, 5)<p t ^ s> < 2 < &- itS ^, H D (e, S)£e 



itH(e,5) 



1p > + 



1 



-UH{e,8) 



+ 2 < ftk 



-Ae-ive-p+e^D(£,5) 



Using again equation and the definition of £, we get 



\<^,~A^>\ < -, 



\<£ip,-iV£-pip >\ < 
\<£iP,ZH D (e,6)iP>\< 



C 



£ l/4j3/2 : 

c 



1p > 



and corresponding equations with e~ ltH ^' 5 ^ip instead of ip. If we suppose that 
the sequence 8(e) satisfies (|2Tj) . then all the terms grow at most as e -1 , so 
we obtain in the end 



< (j) t ,e,S,H D (e,S)(f)t,e,S >< 



c 



Repeating the proof of 



we can then show that 
|| Ce -^M)^ || < CeW^ 



and substituting this back in (jHZJ) we get 



■>t,eM 



< c 



1 1 \ e 1 ' 2 

+ 



5 3/2 £ l/4 ' S 2 5 



c 



.1/4 £ l/2 
+ 



£5/2 £3 



(38) 



□ 



Now, let us fix a sequence 8(e) satisfying (J2HJ). As in last section, we 
consider the time evolution of a product state localized inside one of the two 
edges, away from the vertex, 



-itH(s,5(e)) 



A, 



(39) 
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where (x, y) is the system of coordinates associated to one of the edges, 
/ G C^°(R) and $^ has been defined in (JHJ. If we choose e sufficiently small, 
the tubular coordinates associated to the curve T^s, (ss,u$), coincide with 
(x, y) apart from a small neighbourhood of the vertex. The state ipo is then 
well defined and independent of 5. The limit e — > gives us therefore the 
leading behaviour of an initial state which propagates through a tube which 
curves slowly with respect to the transverse wavelength. 

Equation (|26|) allows us to discard F dg(e)> s( £ )ip £ 5 while (}28|) allows us to 

approximate F ds{e)<s ^ e t with e' itHD ^ e ' s( - £ ^ F ds(e)<s ^ij; . 
We can now prove 

Proposition 2. Let ipl be given by (|39)1 . then, fort G [0, T], 
\\exj>[-UH D (e,S(e)))F dsie)<s{£) i) Q + 

- exp[-itK{5(e)) - itE n /e](f) ■ F dm<s{£) $ £ n (u s{£) )\\ 0, e -> 0, 
where 

K(6(e)) = -ld!-^f . (41) 

Proof. The proof is an application of the fundamental theorem of calculus 
(also called, in this context, Duhamel formula), 

\\exp[-itH D (e,5(e))}F ds{e)<s{£) ^ + 

- exp[-itK(5(e)) - itE n /e}(f) ■ F ds{E)<S{e) <5> £ n (u s{£) )\\ = 

= \\{exp[itH D (e,8(e))]exp[-itK(8(E))-itE n /e]-I}F d5ie ^^ 

+ F dm<s(£) f(x)%(u s ) - F ds(s)<5(£) f(x)<S> £ n (y)\\ < 

< \\{exp[itH D (e, 5(e))} exp[-itK(5(e)) - itEje] - I}F ds(e)<s{e) f ■ $ £ > 5 )|| + 

+ \\F ds(e)<s(£) f(x)^ e n (u 6 ) - F dm<5{£) f(x)$ e n (y)\\ = 

= \\{ew[itH D (e,5(e))]eaq>[-itK(6(e))-itE n /e]-ty^ 

because the second term is zero if e (and therefore 5) is sufficiently small. 
Applying now Duhamel formula 1 , we have that 

\\{exp[itH D (e, 5(e))} exp[-itK(8(e)) - itEje} - l}F ds{e)<s{£) f ■ <S> £ n (u s )\\ < 

<[ ds\\[H D (e,5(e))~K(5(e))-E n /e}- 
Jo 

■ exp[-isK(5(e)) - isE n /e}F ds{e)<&(e) f ■ § £ n (u s )\\. 

(42) 

1 For every fixed <5, the domain of K(8) is H 2 (R), so exp[— itK(6)}(f) ■ F d6< $$(us) is in 
the domain of Hd(s, 6). 
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The Hamiltonian Hjj(e, 6(e)) in curvilinear coordinates, acting on L 2 (K x 
[0, 6), dsdu), is given by 

^ 6{e)) = -\ [i + uk s {s)]* l& + [IT^RF ^ (g) £ + v(s ' u)+ 

- l -^ + ^u 2 

2 5m 2 2e 2 ' 
where V is the geometric potential, 

l ' ' 2\ 4[1 + ufcj] 2 2[1 + uk s f 4 [1 + u* s ] 2 J' 1 ; 

Making a unitary dilation by the factor e 1//2 along u, we get an operator 
acting on L 2 (R x [0, 5(e)/e 1 ^ 2 ], dsdu), given by 

D £l/ ,H D (e,S(e))Dl 1/2 = -\ (1 + J, uh[s)? ^ + |^^T7^^p^ 1/2 ^(s)^+ 

(44) 

where 

Therefore, equation becomes 

||{exp[ittf D ( £) 5(e))] exp[-it£(<f(e)) - itEje] - I}F ds(e)<s{E) f ■ & n (u s )\\ < 

< [ ds \\[D sl/2 H D (e,5(e))Dl 1/2 -K(5(e))-E n /e\- 
Jo 

■ exp[-isK(5(e)) - isE n /e]f ■ F dm<s{e)/E i/2 ■ $* = V«)II ■ 

Therefore, it is clear from previous equations that 

\\[V(s,e l / 2 u) + k 2 /8] eM-isK(6(e)) - isE n /e]f ■ F dm<m/£l/2 • ® E n =1 (u 5 )\\ = 
= 0(e 1/2 /5 3 ) -»• 0, e -> 0, 

' +— u 2 -— )exp[-is^(«J(e))-is£7 n /e]/- 



2e du 2 2e 

■ F d s(e) <S(e)/e 1 / 2 ■ ®T ( U t) 



0, e^0, 
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so we need to control only the terms containing the derivative with respect 
to s in (JED. 

Using lemma El proved below, we have that 

_J_ £ v 2<(s) | exp[ _ i(j?We)) _ itEn/e]f{s) . 



■ F d s(e) <5{e)/eV 2 ■ ®n ( U s) 



< (e^\\k> 



6\\L° 



\dJ\\ L * + \t\-eV^k 5 \\ L oo\\k' s \\lc 



\u&=\ u )\\=0(e 1 ' 2 /5 b )^Q, e^O. 
In the same way we have also that 



2 1 [1 + eV 2 uk s (8) 



ds'- 



exp[-itK(5(e)) - itEje]f{s) 



■ F, 



1/2 • $ e -\u S 



^Ce^HMooIII^/ll + M- 



d s{E) <S(s)/e 

(2 1| A;«5 1| oo • ||^||oo||<9 s /|| + Halloo II /II + Halloo Halloo 



0(e l/2 /5 5 ) -> 0, 
□ 



Remark 4. As stressed above, in this model the dynamics under a strong 
constraining potential is well approximated by Dirichlet boundary conditions 
on a "large" tube surrounding the smooth curve. Proposition El says that 
this choice gives the same result as the procedure of constraining first the 
particle to the motion to the curve, and then taking the limit when the curve 
approaches the graph. 

Lemma 3. Let H be the one- dimensional Hamiltonian H = — ~d x +V, where 
V is a potential bounded together with its first two derivatives, then, given 
ip e H 1 ^), we have 



(45) 



(46) 



\\d x exv(-itH)ip\\ L2m < \\d x ip\\ L 2 m + \t\ ■ \\d x V\\ L *° 
and, given if G iJ 2 (IR) ; 

||^exp(-z^)^|| I/2(R) < ||d^lU 2 (R)+ 

+ l*l(2||9 x ^IU« • ll^vlU 2 w + II^IU« • IMUW 

Proof. Since V is bounded, the domain of the quadratic form associated to 
H is if 1 (IR), and the time evolution sends it into itself. It makes therefore 
sense to write, for ip 6 iJ 1 (IR), 

"t 



[-id x ,e-« H ty = e 
ft 



-itH 



ie 



itH 



ds d s e isH 
ds e lsii [H,-id x ]e- isS tP = e 
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id x )e 



-itH 



ds e isH d x V e~ lsti ^ 



H , 



but this implies immediately 



- id x e~ itH i> = e- itH {-id x )^ + [-id x , e- itH ]iP 

=► ||-i^e-^|| < W-idM + / ds H^Ve-^VH, 

Jo 

which gives ()45|) . 

Following the same path and noticing that 

[H, -d 2 x ] = -d 2 x V - 2d x V ■ d x (47) 

we get dllj). □ 

To complete the analysis of this case we need to study the limit of the 
dynamics exp[— itK(5(e))] when e — > 0. 

The limit of one-dimensional Hamiltonians containing rescaled potentials 
has been studied in detail in the context of the approximation of singular 
interactions, like the delta coupling, by short range smooth potentials (Al- 
beverio et al. (2005) and references therein). The scaling used by us in (}4*Tj) 
however, is not covered in the results presented in Albeverio et al, but it can 
be analyzed using exactly the same techniques. 

The idea is to show convergence in norm of the resolvent of K(5(e)) to the 
resolvent of the Hamiltonian with Dirichlet boundary conditions in s = 0. 
As it is well known (Reed and Simon (1972), theorem VIII. 21) this implies 
strong convergence of the corresponding unitary group. 

One could expect convergence to Dirichlet boundary conditions because 
the potential —kg/8 is a strongly attractive well, which becomes deeper and 
deeper, but whose range is shorter and shorter. As explained in Englisch 
and Seba (1986) in a different context, we expect this to give rise to Dirichlet 
boundary conditions. This in particular says that the strong convergence of 
the unitary group (or the norm resolvent convergence) does not capture the 
behaviour of the eigenvalues which go to — oo when 5 —>■ 0, because, even 
though the ground state of K(5(e)) tends to — oo, its resolvent converges 
to that of a semibounded operator. This phenomenon has already been 
illustrated in Gesztesy (1980). 

We can now prove 

Theorem 4. Let [K(5(e)) — z 2 ]^ 1 be the resolvent of K(5(e)), where ^sz > 0, 
then 

[K(5(e))-z 2 ]- l ^[K D -z 2 y\ e^O, (48) 

in the norm of bounded operators on L 2 (M) ; where Kd is the free Laplacian 
on L 2 (R) with Dirichlet boundary conditions in s = 0. 
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Proof. The potential Q(s) := — k 2 /8 is in L : (R) fl L°°(R), so we can apply 
the dilation technique described in Albeverio et al. (1984) (see also Albeverio 
et al. (2005)). Applying lemma A.l of Albeverio et al. (1984) we get 

[K(5(e)) - z 2 ]- 1 —G z — A s (z)[5 + B s {z)]- 1 C 8 {z) ) Zz > 0, (49) 

where G z is the free resolvent, with kernel g z {w), 

G z :=(H -z 2 )-\ 0»:=^-e fa M 

8 2 2Z (50) 

H := - — , D(H ) = H 2 (R), 

while As(z), Bg(z) and Cs(z) are Hilbert-Schmidt operators with kernels 

A s (z,s,r)=g z (s-5r)\Q{r)\ 1 / 2 , 

B s (z, s,r) = -\Q(s)\ 1/2 9z [5(s - r)]\Q(r)\^, (51) 
C 5 (z,s,r) = -\Q(s)\ 1 / 2 g z (5s-r). 

It is not difficult to see (lemma 2.3 Albeverio et al. (1984)) that 

A s -> Aq, 

Bs -> Bo, (52) 
G$ — > Co, 

in Hilbert-Schmidt norm, where A , S and Co have kernels 
A (z,s,r)=g z (s)\Q(r)\ 1 / 2 , 

B (z,s,r) = -g z (0)\Q(s)\^ 2 \Q(r)\ 1 / 2 , (53) 
C (z,s,r) = -\Q(s)\^ 2 g z (-r). 

The operator B is not invertible on the whole Hilbert space, but it is clear 
from the expression of the kernel that it actually acts on the one-dimensional 
subspace, denoted by Hq, generated by the vector tpq given by 

So we can write 

B = -gz(0)\\Q(s)\\ L i m (p Q <<p Q ,-> . (55) 
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On Hq, B is invertible and the inverse is given by 



B o 1 = fn\\\nf mi < ¥>Q» • > ■ (56) 

9z{0)\\Q{s)\\ L i m 

Since the operator C has range equal to Hq and A acts non trivially 
only on Hq, we get that 

[K(S(e)) - z 2 }' 1 - G z - AoB^Cq, 
which has a kernel given by 

a ( s r ) 9z(s)g z (-r) 

g z {s r) gM , (57) 

which is the kernel of the resolvent of the Dirichlet Hamiltonian. □ 
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